Spectroscopic and thermodynamic properties in a four-band model for pnictides 
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In this paper we provide a comprehesive analysis of different properties of pnictides both in the 
normal and superconducting state, with a particular focus on the optimally-doped Bai_a;Ka;Fe2As2 
system. We show that, by using the band dispersions experimentally measured by ARPES, a four- 
band Eliashberg model in the intermediate-coupling regime can account for both the measured 
hierarchy of the gaps and for several spectroscopic and thermodynamic signatures of low-energy 
renormalization. These include the kinks in the band dispersion and the effective masses determined 
via specific-heat and superfluid-density measurements. We also show that, although an intermediate- 
coupling Eliashberg approach is needed to account for the magnitude of the gaps, the temperature 
behavior of the thermodynamic quantities does not show in this regime a significant deviation with 
respect to weak-coupling BCS calculations. This can explain the apparent success of two-band BCS 
fits of experimental data reported often in the literature. 

PACS numbers: 74.20.-z, 74.25. Jb, 74.25. Bt 



I. INTRODUCTION 

The recent discovery of superconductivity in pnictides^ 
prompted an intense experimental and theoretical re- 
search about the properties of these materials. At 
the very beginning the analogies between pnictides and 
cuprate superconductors (e.g. the layered structure and 
the phase diagram) suggested that a similar route to 
high-temperature pairing could be at play in these two 
classes of materials. ^''^ However, a large experimental ev- 
idence has been accumulated so far that significant dif- 
ferences between pnictides and cuprates are also impor- 
tant, starting from the very basic fact that pnictides have 
a multiband structure. According to LDA calculations, 
indeed, the band structure of pnictides near the Fermi 
level is characterized by two hole-like bands around the 
r point, and two electron-like bands around the M points 
of the reduced Brillouin zone.'*"'' A third hole- like band at 
the r point could be expected to cross the Fermi level in 
some materials, but eventually it moves below the Fermi 
level when the experimental value of the apical As posi- 
tion is used in LDA calculations.^ 

Despite the large theoretical work devoted to address 
the outcomes of multiband superconductivity, many open 
issues still remain about a direct comparison between 
the theoretical predictions and the experiments, or be- 
tween the outcomes of different experimental probes. A 
first issue concerns the experimental observation of only 
two gap values in hole-doped 122 compounds,*'^ whereas 
in a multiband BCS approach one would generically ex- 
pect a different gap in each band, depending on the cou- 
pling and on the density of states (DOS) of the sev- 
eral pockets involved in the pairing. ^''""'^^ This is true in 
particular for hole-doped Bai„xKa;Fe2As2(BKFA), where 
many detailed experimental findings have been accumu- 
lating due to the existence of large crystals. Here angle- 
resolved photo-emission spectroscopy (ARPES) has re- 



ported quite different DOS in the hole and electroir pock- 
ets involved iir the pairing,^ so that at the BCS level one 
could expect to observe three gaps, one for each hole band 
and one for the (almost degenerate) electron bands. On 
the contrary. ARPES experiments have reported only two 
different gaps: a large one on the inner hole pocket and 
on the electron ones (with A/Tc ^ 3.5), and a small one 
(with A/Tc ^ 1.8) on the outer hole pocket. ^'^ 

Besides the iroir-BCS hierarchy of the gaps, further 
difficulties arise in the attempt to reconcile ARPES 
data with several thermodynamic measurements. For in- 
stance, photoemission experiments performed by several 
groups in different pnictide materials have shown that 
there is a substantial renormalization of the whole band 
structure with respect to LDA predictions, with a reduc- 
tion at least of a factor two.^'^^"*^ At the same time, the 
estimates of the specific heat coefficient Cv/T obtained 
by using the ARPES bandwidth, despite being substan- 
tially larger than LDA, are still about a factor two smaller 
than the values measured in the normal state for 122 
compounds. ^^"^® This comparison calls for a dichotomy 
between high-energy and low-energy mass renormaliza- 
tion, that must be accouirted for by different mecha- 
nisms. Recently, a similar distinction between renor- 
malization effects operative at different energy scales has 
been pointed out also in optical-conductivity measure- 
ments in 1111 compounds.*^ 

As far as the temperature dependence of the specific 
heat in the superconducting state is concerned, the com- 
parison with ARPES is again compelling: indeed, de- 
spite the large A/Tc values reported by ARPES that 
call for an intermediate/strong coupling pairing mech- 
anism, the temperature profile of the specific heat can 
be well reproduced by a simple BCS fit.'^^'^^ A simi- 
lar result arises from the airalysis of superfluid-density 
measurements, ^°"^^ where two-gaps BCS fits seem to 
work quite well once that the experimental A/Tc ratios 
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are implemented. 

In this paper we provide a systematic analysis of spec- 
tral and thermodynamic properties in pnictides with the 
goal of reconciling the results obtained with the differ- 
ent probes. Our analysis is based on a four-band model 
where carriers interact with bosonic excitations treated 
within the Eliashberg approach. As we shall sec, the ob- 
served hierarchy of the gaps calls for a predominant inter- 
band nature of the interactions, making spin fluctuations 
the most natural candidates for the pairing glue.^'''^^ We 
shall focus in particular on the effects of the exchange 
of spin fluctuations on several spectroscopic and thermo- 
dynamical properties. Since the typical energy scale for 
spin fluctuations is of the order of 20 meV,^^~^^ it cannot 
be responsible for the overall band narrowing observed 
by ARPES, that is operative up to a rather high energy. 
In this context we shall extract input band parameters 
directly from the high-energy ARPES measurements, by 
assuming that their renormalization with respect to LDA 
calculations originates from electronic correlations. 
This approach is thus different from what discussed in 
Ref. [11], when the experimental ARPES determination 
of the bands for doped BKFA was not yet available. In 
that case it was shown that, within a BCS approach, the 
experimental observation of only two gap values could be 
accounted for by a suitable (moderate) renormalization 
of the LDA band parameters. Here instead we show that 
the experimental measurement of the band structure to- 
gether with the gap values on each band provide in this 
material a compelling constraint for the microscopic the- 
ory. Within this framework we estimate the magnitude 
of the different interband couplings from a comparison 
with the measured gaps. We find that the dimensionless 
couplings vary from A~0.2toA~ 1.6, depending on the 
band. We also calculate the additional mass renormal- 
ization due to the exchange of spin fiuctuations. These 
low-energy features are hardly visible in ARPES but they 
are responsible for the large effective mass of the charge 
carriers probed by specific-heat measurements, that are 
sensible to excitations near the Fermi level, solving then 
the apparent contradiction between the different set of 
measurements. Finally, we analyze the temperature de- 
pendence of the specific heat and superfiuid density, and 
we show that at the coupling values relevant for pnictides 
we do not observe significant deviations fi-om a conven- 
tional BCS profile, explaining the apparent success of the 
BCS fits proposed in the literature. 

It is worth pointing out the differences between our 
approach and previous works on multiband Eliashberg 
calculations proposed in the literature. ^^'^^ The tendency 
of the gaps to assume the same value in strongly nested 
bands within the Eliashberg theory was already noted in 
Ref. [12]. However, the authors considered there a two- 
band model, so that it was impossible to reproduce the 
second smaller gap value measured by ARPES, which is 
realized in a third, less coupled band. Indeed, a cor- 
rect approach to pnictides requires using at least a four- 
band model with an anisotropic interband pairing, as 



it was pointed out previously within a BCS scheme in 
Ref. [11]. An Eliashberg approach to a four-band model 
has been explored recently in Ref. [36], where the au- 
thors were aimed to reproduce exactly the experimental 
ratios A/Tc in the various bands. An extremely large 
coupling A > 4 was there found for BKFA. Such anal- 
ysis disregards however the fact that an accurate esti- 
mate of Tc within the mean-field-like Eliashberg theory 
is doubtful in these almost two-dimensional materials, 
where superconducting fiuctuations are expected to be 
relevant due to the low dimensionality,^'^'^^ leading to a 
lowering of the real T,, in comparison with the mean- 
field estimate. In this situation we prefer to concentrate 
our analysis on the consistency between the gap values 
and the density of states. As mentioned above, assum- 
ing a typical energy scale ~ 20 meV for the charac- 
teristic spin-fiuctuations, we get Ai ~ 0.2 — 1.6. much 
lower than in Ref. [36] (i being here the band index). 
On the other hand, as we shall show below, these values 
appear to be perfectly compatible with the thermody- 
namical properties, whereas stronger couplings would be 
inappropriate, because the low-energy renormalization of 
the charge carriers would be too large compared to the 
experimental outcomes from specific-heat and superfluid- 
density measurements. Our estimates of the interband 
coupling Aj ~ 0.2 — 1.6 in BKFA locate this material in 
the weak-intermediate coupling regime. This observation 
could suggest that analytical expressions d la McMillan- 
AUen-Dynes^^ would be appropriate, as proposed in Ref. 
[12]. This is however not the case in multiband systems 
where, as we show below, McMillan-Allen-Dynes-like ex- 
pressions can qualitative fail already above very weak 
coupling A > 0.2, so that a numerical solution of the 
multiband Eliashberg equation is required. 

The structure of the paper is the following. In Sec- 
tion II, we briefiy review the results of a two-band model 
in order to elucidate the differences between the BCS 
and Eliashberg approach and the need of a numerical 
solution even in the weak/intermediate coupling regime. 
The reader interested only in the comparison with the 
experiments can skip this technical discussion, and refer 
directly to Section III, where we introduce the full four- 
bands model, and we show that at intermediate coupling 
the Eliashberg theory can reproduce the experimentally 
measured gap values in pnictides. In Section IV we show 
the results for the specific heat and the superfiuid den- 
sity. Finally, in Section V we draw some conclusions and 
we discuss the perspectives of our work. 



II. TWO-BAND MODEL 

The BCS theory is characterized by a number of uni- 
versal behaviors (as the Tc vs. A relation, the A/Tc — 
1.76 ratio, etc.) which are strictly valid only in the 
limit where the dimensionless coupling A — )• 0. Devia- 
tions from these universal results are related to interme- 
diate/strong coupling effects, so that the analysis of such 
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deviations could be employed in principle to estimate the 
strength of the coupling A. McMillan-AUen-Dynes-like 
formulas,^® based on a controlled expansion in power of 
A, could be quite useful in this context because in single- 
band models they provide analytical expressions to quan- 
tify these effects without resorting to a numerical solution 
of the Eliashberg equations. 

In this Section we show that in the multiband case with 
predominant intcrband interaction the situation is quite 
different. Indeed, the McMillan-Allen-Dynes-like expan- 
sion reproduces the Eliashberg behavior as function of A 
only for very weak coupling A < 0.2, whereas for larger 
couplings a full numerical solution of the Eliashberg equa- 
tions is required. We demonstrate this result for simplic- 
ity within a two-band system, previously addressed in 
Ref. [12]. We thus write the general Eliashberg equa- 
tions for a purely interband interaction which is taken to 
be repulsive in the Cooper channel: 



Zi{n)Ai{n) = -Ai27rT^ Z)(n-m) 



A2(m) 



(1) 



Z,{n) = l + Ai2^^Z)(n-m)^== 



(2) 

together with an equivalent set of equations for A2 , Z-2 
related through a coupling constant A21. Here for sake 
of shortness we denote the dependence on the Matsubara 
frequency w„ of the gap function A^ and the the renor- 
malization function Zi for the band i as Aj(n) = Aj(ia;„) 
and Zi{n) = Zi{iujn). D{n ~ m) ~ D{ujn — uJ„i) is 
the boson propagator, which is related to the Eliash- 
berg spectral function i?(r2) by the relation D{n — m) = 
/ 2rtdnB{n) /[{ujn - ujrn)'^ + n'^]. The dimcnsionlcss cou- 
pling constants A12, A21 can be expressed in term of an 
unique energy coupling G > weighted by the appropri- 
ate density of states A^,, namely A12 = GA''2, A21 = GA^i. 

Eqs. (l)-(2) can be solved self-consistently to obtain 
a numerical exact solution of the Eliashberg equations, 
assuming, for simplicity, an Einstein boson spectrum 
B{fl) = {uJo/^)S{0, — Wo), where ojq is the characteris- 
tic boson energy. For a repulsive interaction, G > 0, the 
gaps in the two bands have opposite signs, so that the or- 
der parameter has a s± symmetry. In the rest of the Sec- 
tion we will assume conventionally Ai > and A2 < 0. 
Moreover, to make a direct comparison with Ref. [12], we 
consider the case where the ratio B = N2/N1 of the DOS 
in the two bands is B = 2.6. 

Let us focus first on the gap anisotropy A = A1/IA2I 
at T = as a function of the average coupling A = 
■V/A12A21, which was extensively analyzed in Ref. [12]. 
As one can see in Fig. la, A \fB as A — >■ 0,^^ 
but within the Eliashberg framework A approaches 1 as 
A increases, showing that the gaps get closer to each 
other. This result is in sharp contrast with the BCS 
solution, that is obtained from Eqs. (l)-(2) by ne- 
glecting the equation for the renormalization functions 



\Zi(ri) = 1], and assuming a BCS factorized square- well 
model D{n ~ m) = 9(ujq — |a;„|)0(u;o — Within 
this framework A,(n) = Ai^(a;o — \^n\) and one gets the 
simple equations 



Ai = -Ai2A2n2, 

A2 = -A2iAini, 



(3) 
(4) 



where = 7rT^„0(wo - |w„|)/ y^of+Af. The be- 
havior of A obtained by the numerical solution of the 
previous BCS set of equations is also reported in Fig. la: 
as one can see, the two gap values diverge one from the 
other as the coupling increases, in contrast to the results 
of the intermediate-strong coupling Eliashberg solutions 
of Eqs. (l)-(2).i2 

Since the Eliashberg theory accounts for the effects of 
the Z-renormalization functions, in Ref. [12] it was pro- 
posed a simple analytical way to illustrate the difference 
between BCS and Eliashberg approach by means of a 
"renormalized BCS model". In this case the square-well 
model for the gap equations can be completed with a 
corresponding square- well model for the renormalization 
spectral functions, Zx(n) = 1 + A12, Z2{n) = 1 -h A21 for 
|w„| < cjo and Zi{n) = Z^in) = 1 for > uq, so that 



AiZi 
Zi 

A2Z2 
Z2 



— Ai2A2n2, 

1 + A12, 
-A2iAini, 
1 + A21. 



(5) 
(6) 
(7) 
(8) 



The gap anisotropy obtained from the above set of equa- 
tions is also shown in Fig. la, compared to the Eliash- 
berg and BCS solutions. Remarkably, one sees that 
the renormalization effects account very well for the de- 
creases of the gap anisotropy A as the coupling constant 
increases, giving essentially the same A{X) dependence 
as the Eliashberg calculations. As it was suggested in 
Ref. [12], this result can be understood by an analytical 
approximation of equations (5)- (8) at low coupling. In- 
deed, at T = one can use the approximate BCS forms of 
the Ui bubbles, Ui = sinh"^ (wo/jAil) w log (2a;o/|Ai|) 
to write a self-consistent expression for the ratio A as a 
function of the dimensionless couplings A12 and A21: 



AZi 



12 



AX 



21 



log^. 



(9) 



As a consequence, Eq. (9) can be solved perturbatively 

in powers of the effective coupling A = Xj^fZiZ^-, which 
takes into account, at a BCS level, the self-energy renor- 
malization: 



(10) 



A = VB{1 + cX + dX^) + 0{y), 



where B = {Z2/Zi)B. By substituting Eq. (10) into Eq. 

(9), and recalling that Zx\f^lXx2 = Z2/\^X2i = 1/A, 
we immediately obtain: 



c=-logB, d=^— 



41ogB-|-log^B 
32 ■ 



(11) 
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While making an expansion in the coupUng constant one 
must properly expand also A and B at the leading order 
in A, namely, 



B(l + A2i -A12) 



B 
A 

As a consequence, one finds^^ for the gaps ratio (10) 



(12) 
(13) 



A 



Ai 
IA2I 



Vb 1 



A_ 
~2~ 



A 



■logB 



(14) 



where A± = A21 ± A12. The BCS case is recovered from 
the previous equations by setting Zi = Z2 = 1, i.e. 
X± = 0. Thus, since A_ < when B > 1, and the 
second term is larger than the third one regardless the 
value of B > 1, from Eq. (14) it follows that in the BCS 
case A increases with increasing A, while in the presence 
of renormalization effects A decreases, i.e. the two gap 
values approach each other, as confirmed by the numeri- 
cal solutions at all A values reported in Fig. la. 

From the above considerations and the results of Fig. 
la one could then be tempted to conclude, as it was 
done in Ref. [12], that the McMillan- Allen-Dynes-like 
equations (5)-(8) capture the basic physics of the Eliash- 
berg solution at all coupling values. However, this is 
not the case, as we show in Fig. lb, where we report 
explicitly the Aj/Tc values in the various approaches. 
Here we use for simplicity in the BCS and renormalized- 
BCS case the estimate Tc = 1.13wo exp (— 1/A) and 

Tc = l.lBwoexp ^— 1/A^, respectively, valid at weak cou- 
pling by means of the approximate BCS form H^ = 
log (l.lBwo/Tc) of the bubbles near Tc. As one can see in 
Fig. 1, in the renormalized BCS case the gap values ap- 
proach each other by a decrease of the larger Ai /Tc value, 
and a partial increase of the smaller |A2|/Tc value. This 
result can be again understood analytically at low cou- 
pling by resorting to the above expansion (10) and the 
Tc expression. We then obtain the leading dependence of 
Ai/Tc on the coupling: 



Ai 

Tc 



1.76i?V^fl + A^^°g^-^"g^^ 



Tc 



32 



™^/^(i + ^ + All^M!^)a5) 

1.76i^-V^(l-A^^°g^3+^"g^^) 



l.,g^-V4^1_A^_,41og5 + lo^ 
V 4 32 



16) 



Also in this case the terms in A_ < are larger than the 
others, so that one recovers from the above equations that 
at low coupling Ai/Tc decreases and |A2|/Tc increases as 
a function of A. However, except for a very narrow range 
of coupling A < 0.2 — 0.3, the numerical Eliashberg solu- 
tion of Eqs. (l)-(2) is markedly different. Indeed, in the 
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FIG. 1: (Color online) Evolution of the ratio A1/IA2I (panel 
a) and of |Ai|/rc (panel b) as function of A calculated by solv- 
ing numerically the BCS Eqs. (3)-(4), the renormalized BCS 
Eqs. (5)-(8) and the Eliashberg Eqs. (l)-(2), with B = 2.6. 
Notice that while the renormalized BCS model reproduces 
the behavior of A1/IA2I in Eliashberg theory at all coupling 
values, the same is not true for the |Ai|/Tc values. 



intermediate/strong coupling regime, which is the one 
relevant for pnictides, both Ai /Tc and |A2|/Tc increase 
with A in the Eliashberg case. Thus, in the full numerical 
Eliashberg solution the gaps in the two bands approach 
each-other by means of an increase of the absolute Ai/Tc 
ratio in both the bands, that is the typical signature of 
strong coupling. We note in passing that while Fig. la 
reproduces the findings of Ref. [12], the same is not true 
for Fig. lb. Indeed, the numerical results for Ai/Tc in 
the Eliashberg theory reported in Ref. [12] differ signifi- 
cantly from our findings, even in the low-coupling regime 
where the above analytical analysis supports completely 
our numerical calculations. 



III. SUPERCONDUCTING PROPERTIES IN A 
FOUR-BAND MODEL 

In the previous section we have shown within a sim- 
ple two-band model that a full numerical approach is 
needed to capture the property of Eliashberg equations 
of removing the DOS anisotropy of the gaps in the pres- 
ence of interband pairing. Following the same reasoning 
we focus in this section on a four-band model, to cor- 
rectly capture the physics of pnictides. In particular, 
we shall discuss the case of hole-doped Bao.6Ko.4Fe2As2, 
using the notation of Refs. [8,9], where a (/3) is the in- 
ner (outer) hole-pocket centered around the F point, and 
7i , 72 are the two electron-like pockets centered around 
the M points of the folded Brillouin zone of the FeAs 
planes (see Fig. 2). The dominant interactions in pnic- 
tides are thought to be mainly interband ones, connect- 
ing hole Fermi-sheets with electron Fermi-sheets, through 
the exchange of spin fluctuations at the antiferromag- 
netic wave-vector Q = (tt, tt).^^'^^ The strength of such 
interband coupling between hole and electron pockets. 
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FIG. 2: (Color online) Schematic view of the four-band model 
we use for pnictides. a and P are the hole bands around the 
r point, 7i — 72 are the two degenerate electron bands at 
the M point. The sizes of the pockets are inferred from the 
experiments in Bao.6Ko.4Fe2As2.^ The coupling anisotropy, 
with G > g is suggested by the different nesting properties of 
electron to the hole bands, as due to their different sizes. 



in these materials, is in addition expected to depend on 
the relative sizes of the pockets, that naturally affect the 
nesting condition for the spin fluctuations mediating the 
pairing. In particular, in BKFA the size of the a and 
7i Fermi surfaces is quite comparable, while the /3 band 
has a Fermi surface substantially larger, with a corre- 
sponding less degree of nesting. In this situation the in- 
terband /3 — ji coupling g is expected to be significantly 
smaller than the a — 7^ coupling G. In addition, since the 
two electron pockets have comparable sizes, we assume 
them for simplicity to be degenerate, and we denote with 



N. 



N. 



71 



the total DOS in the electron pockets. 
We can write thus the multiband Eliashberg equations in 
the form: 



Za{n)Aa{n) 
Zp{n)Ap{n) 

Z^{n)A.y{n) 



GNaAa{m) 

V< + A2(m) 

(19) 



— ttT^^ D{n — m) 

rn 

^< + A2(m)_ 



IH 'S^ D(n ~ m) 

m 

TtT ^ — ^ 

1 H 2_^D{n ~ to) 



iH y D{n - to) 



__gNpUJrn__ 

^L^l + Alim) 



GN^UJ.m 



^C^2^+A2(to) 
9_NjUJrn__ 
UA2(to) 

Vc.2^ + A2(to) 



(20) 
(21) 



(22) 



In the BCS limit, the above equations reduce to the one 
already discussed in Ref. [11]: 



A„ 
A. 



-N^gA^U^ 

GN^A^U^+gNpApUp, 



(23) 
(24) 
(25) 



and Tc is given by Tc = 1.13woe~^/^, where, in anal- 
ogy with the two-band case, we introduce the effective 
coupling 



A 



^G^N^N^+g^N^Np. 



(26) 



It is interesting to notice that in the BCS limit the 
ratio between the two hole gaps is simply given by 
the ratio of the corresponding interband couplings, i.e. 
Aa{n) / Ap{n) — g/G^ mdepenrfent/?/ on the relative DOS. 
The simple experimental observation Aq,/A^ « 1/2 
would suggest thus, in BCS, g k, G/2. Such con- 
straint does not apply however to a more accurate Eliash- 
berg analysis [Eqs. (17)-(22)], where AQ(n)/A^(n) = 
{g/G)Zp{n)/Za{'n), so that, in principle, the ratio of the 
gaps in the two hole pockets depends both on the cou- 
plings and on the DOS of the various bands, requiring 
thus a more careful analysis. 

We employ now Eqs. (17)-(22) to evaluate the micro- 
scopic interband couplings G, g from the physical con- 
straints given by the experimental determination of the 
gap magnitudes on the different bands. We shall apply 
later this analysis to calculate different superconducting 
and normal-state properties, as the superfluid density 
and the specific heat, in order to have an independent 
check of the reliability of our analysis. 

A debated issue in this context is the assess of a proper 
choice for the underlying normal state electronic bands. 
Indeed, as we mentioned in the introduction, ARPES 
measurements in several pnictide families report signif- 
icant differences in the electronic dispersion compared 
with LDA calculations, with an apparent renormaliza- 
tion of the whole band structure by a factor 2.^'^'^"^^ 
Most striking, such band narrowing seems to be oper- 
ative up to very high energy scales, as it is confirmed 
also by recent optical sum-rule analysis performed in 
LaFePO samples. This overall renornialization of the 
bands with respect to LDA seems thus a general feature 
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e° (meV) 


28 


43 


-60 


ti (meV) 


54 


27.5 


160 


Ni (eV-i) 


1.47 


2.89 


0.50 


Ui 


0.082 


0.24 


0.06 



TABLE I: Microscopic band parameters extracted from Ref. 
[9] by approximating each band with a parabolic form ei{\f.) = 
el — ti|kp. Ni = l/Anti is the DOS and is the number 
of holes/electrons per unit cell in each band. Note that for 
the 7i , 72 bands we accounted for the corrections due to the 
elhptical shape reported in Ref. [9]. 



of pnictides, probably arising from local Hubbard-like 
correlations, '^^"^^ so that it cannot be captured by the 
coupling of the electrons to low-energy bosonic modes. 
To take into account this feature, we estimate our in- 
put band parameters directly from ARPES experiments, 
that have enough resolution to capture the high-energy 
mass renormalization. Using the tight-binding fit of the 
bands reported in Ref. [9], we approximate close to the 
Fermi level each band as eiCk) = — ti|kp, where k is 
measured with respect to the F point for the hole bands 
and to the M point for the electron bands (see Fig. 2). 
Ni = l/Anti is the DOS (per spin) in each band. The 
band parameters for each band, extracted from Ref. [9], 
are listed in Table L To model the spin-mediate inter- 
action, following Ref. [40] and the recent report [32] we 
use B{U) = ^}ujo/tt{ujq + fi^), with the characteristic en- 
ergy scale ujQ — 20 meV estimated from experimental 
measurements. ■^^"^^ 

In approaching a numerical analysis, we should note 
that, neglecting for the moment the weaker coupling 
7-/3, the model is equivalent to two-band case discussed 
in the previous Section, with B — Na/N^ ~ 1.5, so 
that at the BCS level the gap values in the two strongly- 
nested bands are different at least by 20%. The presence 
of a finite scattering /3 7 makes this scenario even 
more complex, with the onset of the gap A^, which con- 
tributes to increase the magnitude of A-^ . In the absence 
of any particular degeneracy between the a and 7 bands, 
this situation would result, within a weak-coupling BCS 
framework, in the appearance of three different gaps on 
the different Fermi sheets, as marked by the dots in Fig. 
3. The experimental observation from ARPES of two 
nearly degenerate gaps on the a and 7 bands calls then 
for further investigation. 

In this regard, the strong coupling results from the 
two-band model, discussed in Ref. [12] and readdressed 
more specifically in Section II, shed a new interesting 
light once plugged in a four-band model. Within this 
context, indeed, the nesting-driven strong coupling in- 
teraction between the a and 7 bands leads to a merging 
of the value of the two large gaps Aq,, A-^ in the systems. 
At the same time, the anisotropy of the coupling between 
the two hole pockets and the electron one, guaranteed by 
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FIG. 3: (Color online) Dependence of the ratios Aa/Tc, 
Ap/Tc and A^/Tc as functions of the coupling G (g — G/2) 
within the Eliashberg theory. The dots at G = mark the 
weak-coupling BCS values. 



the different nesting conditions, allows the system to keep 
the gap in the /3 band smaller even within the Eliashberg 
approach, in agreement with ARPES. 

To have a more quantitative insight we plot in Fig. 
3 the gap values A^ — Ai{n = 0) obtained from the 
numerical solution at T = Tc/SO of Eqs. (17)-(22) as a 
function of the coupling G, for the indicative case g — 
G/2, which would gives the experimental value A^ = 
0.5 Aq in the BCS limit G — 0. As we can see, the use of a 
four-band model is crucial to recover the hierarchy of the 
gaps: while in the two-band case by increasing A one is 
forced to have a single gap value, in the four-band model 
the anisotropy of the couplings, which follows naturally 
from the different nesting properties of the various bands, 
allows for a two-gap result, where two almost degenerate 
larger gaps are predicted on the electron and on the inner 
hole pocket, while a smaller gap is found in the outer hole 
pocket. 

As mentioned in the introduction, whereas the exper- 
imental gaps in the a and 7 bands are approximately 
equal, the ratio Aa/Tc « 3.75 is quite larger than the 
BCS value 1.76. Recently, the possibility to reproduce 
exactly these ratios in BKFA has been investigated in 
Ref. [36], where it was shown that an extremely large 
value of the effective coupling A > 4 and a small bo- 
son energy scale coq = 10 meV were required. This kind 
of analysis needs however to be taken with some cau- 
tion. Indeed, any estimate of Tc in a mean-field-like 
theory as the Eliashberg one is in general questionable 
in two-dimensional systems, as pnictides, where the su- 
perconducting fluctuations, when properly taken into ac- 
count, could significantly reduce Tc.^^ For this reason, 
while the low-temperature gap values estimated within 
a mean- field- like Eliashberg theory can be quantitatively 
sound, the estimate of Tc done within the same approach 
must be taken as an upper limit. Moreover, the values 
of LdQ and of the coupling do not influence only the ra- 
tios Ai/Tc, but they also control in a crucial way other 
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a /3 71.72 



A, (meV) 


9.48 


4.35 


-10.48 


Z, (meV) 


2.09 


1.35 


3.67 


m* / me 


9.61 


12.28 


5.72 


JsAO) (meV) 


4.8 


10.7 


6.14 


AfP (meV) 


12 ± 1 


6± 1 


12 ± 1 



TABLE II: Eliashberg parameters obtained by numerical so- 
lutions of the Eqs. (17)-(22) with the coupling matrix (27). In 
the last row we report for comparison also the experimental 
values of the gaps from Ref. [8,9]. 

physical quantities, like the mass renormalization and the 
position of the kink in the band dispersion, that can be 
accessed experimentally. For these reasons we investigate 
here the following approach: we discard the exact deter- 
mination of Tc and we choose the coupling strength as 
to reproduce the magnitude of the different gap values. 
Afterwards, we check if by means of the same parameter 
values we can account for other experimental results re- 
lated to mass-renormalization effect in the spectral and 
thermodynamic properties. 

In Table II we summarize our results for the interband 
scattering, G = 1.1 eV , g = 0.32G' ~ 0.35 eV and luq = 20 
meV, obtained to reproduce the experimental gaps. The 
multiband matrix of the coupling constants, in the band 
space (a, /3, 71 , 72) reads: 

/ 0.55 0.55 \ 

A= 0.18 0.18 

1.62 1.01 

\ 1.62 1.01 / 

and the effective average coupling defined in Eq. (26) 
is A = 1.5. We notice that a better agreement with the 
experimental gap values could be enforced by slightly dif- 
ferent DOS than the ones estimated in Ref. [9] , or by as- 
suming a non-zero interband coupling, as due to phonons. 
However, we prefer here to use a minimal set of free pa- 
rameters to show the overall quantitative agreement be- 
tween our approach and the experiments. We obtain a 
critical temperature Tc = 0.21cl'o = 48.7 K, which over- 
estimates the experimental one T^^p = 37 K by about 
10 K. Taking into account that a recent analysis of para- 
conductivity has shown that this is the typical range of 
temperature where superconducting fluctuations are ac- 
tive in pnictides,'^* one can expect that the effect of su- 
perconducting fluctuations beyond Eliashberg theory will 
contribute to improve the agreement between the present 
estimate of Tc and the experimental value. 

The coupling of carriers to spin fluctuations reflects on 
the one-particle spectral properties already above Tc. In 
particular, the bands in the normal state are expected to 
display a kink at an energy wq, so that for E < loq the 
Fermi velocity is renormalized with respect to the bare 
value, with Wp ^ = v-p^i/Zi (or equivalently m* = Zi-m^^j), 




0.1 0.2 0.3 0.4 0.5 0.6 1.3 1.4 1.5 1.6 

k(A-^) k(A-i) 



FIG. 4: (Color online) Intensity map of the spectral function 
in the normal state for the renormalized bands within the 
Eliashberg theory, as compared to the bare band dispersions, 
represented by the dashed lines. Note that within a spin- 
fluctuation model for the bosonic excitation the kink features 
at i5 = —LJo are considerably smeared out in the spectral 
functions. Here we added a constant broadening 77 = 0.002tJo 
to account for a small amount of disorder. 



where Zi = Zi(n = 0) and TOb,i is the band mass of 
the sheet i. In Fig. 4 we show the intensity map of 
the spectral function for the interacting hole and elec- 
tron bands, as obtained by the Marsiglio-Schossmann- 
Carbotte analytical-continuation procedure, along with 
the bare band dispersions. For the spin-fluctuation model 
used here, and for the coupling values deduced by the 
measured gaps, the kinks at ujq are signiflcantly smeared 
out in the spectral function. Unfortunately, the exper- 
imental resolution of the data in Ref. [9] is not high 
enough to resolve the effect of low-energy spin fluctua- 
tions from the high-energy renormalization. However, it 
is interesting to notice that a similar kink has been ac- 
tually observed in high-resolution ARPES measurements 
performed by an other group in a BFKA sample with 
lower doping than the one we are discussing here.'*^ In 
particular, the authors flnd a kink around approximately 
30 meV, and a velocity renormalization in the inner and 
outer hole pockets of about 1.8 and 1.6. Thus, given 
the difference in doping and the lack in our approach of 
a specific momentum dependence of the spin-exchanged 
fluctuations (that can contribute to slightly increase the 
effective energy of the kink), our results are in good agree- 
ment with these flndings. 

Below the superconducting transition the position of 
the kink in the band dispersion is in principle shifted 
in each band to an energy given by ujq + Aj , where Aj 
is the value of the gap in the bands coupled to the i- 
th sheet. In particular, in the a, /3 bands the kink due 
to spin fluctuations is expected to be recovered at an 
energy cjq + Ay « 30 meV, while in the 71 , 72 bands two 
kinks will appear, at the energies wq + Aq « 30 meV and 
cjQ + Ah « 25 meV. It should be pointed out, however, 
that, because of the smearing of the spectra in Fig. 4, 
it would be quite hard to detect the possible splitting 
in the 7 bands due to two different gap values of the 
a and (3 bands. At the best of our knowledge, a clear 



8 



identification below of a continuous shift at higher 
energy of the kinks observed in the normal state has not 
been reported yet in pnictides.^^ 

Finally, it is worth mentioning that an additional 
source of discrepancy between the LDA and the exper- 
imental bands in pnictides comes from finite-band ef- 
fects, that have been discussed in Ref [44]. Indeed, when 
the strong particle-hole asymmetry of pnictides is taken 
into account by considering the finite bandwidth, the 
spin-mediated interband coupling leads to a shift of the 
Fermi momenta with respect to LDA, that has been in- 
deed measured in other pnictides by de Haas- van Alphen 
experiments. ^^'^^ In the present case we did not com- 
pute explicitly these shifts, since they are already in- 
cluded in the band dispersion extracted from the mea- 
sured ARPES data. Since finite-band effects do not alter 
qualitatively the self-energy corrections apart from the 
mentioned energy shift, we solved Eqs. (17)-(22) within 
the usual infinite-band approximation. 



IV. THERMODYNAMIC PROPERTIES 

Having determined an appropriate set of band param- 
eters and multiband couplings, we investigate now the ef- 
fects of the spin-mediated interactions on the thermody- 
namical properties. Indeed, the signatures of low-energy 
renormalization are much more easily detectable in ther- 
modynamic measurements of masses enhancement than 
in photoemission, where a very high resolution is required 
to resolve the kinks in the band dispersion. 

Let us consider as a first insight specific-heat measure- 
ments. In the normal state, the coefficient = Cy jT 
of the linear T term in the specific heat measures the 
mass enhancement, once compared with the value esti- 
mated for non-interacting bands. To clarify the units, 
we shall refer in the following to the specific heat per for- 
mula unit (so that one mole of the materials contains 2 Fe 
atoms). By expressing the renormalized DOS (per spin) 
of each band as a function of the renormalized mass to* , 
N* = TO* /27r and, restoring explicitly all the needed di- 
mensional constants, we have that each band contributes 
to as: 



27r^fc 



2l2 



In.;. 



'•-N*NAa^ = 1.5^ mJ/KVol, (28) 



where a = 3.9 A is the lattice spacing in BKFA, Na 
is the Avogadro number, ks the Boltzmann constant 
and TOg the free electron mass. Within LDA one ob- 
tains = 9.26 mJ/K^mol,^ that is remarkably smaller 
than the values obtained in doped BKFA, either by di- 
rect analysis of the normal-state specific heat 7jv ~ 49 
mJ/K^mol^® or by measurements of the upper critical 
field « 63 mJ/K^mol^^. By means of the band pa- 
rameters extracted from ARPES and listed in Table I 
one estimates 7jv = 25 mJ/K^mol, while with the renor- 
malized masses listed in Table II, that include low-energy 
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FIG. 5: (Color online) Temperature dependence of the differ- 
ence between the superconducting and normal-state specific- 
heat coefficient 7s— 7]v = ACv^/T in the various bands, along 
with the total one. 



renormalization effects on the ARPES bands, we can es- 
timate 7JV — 50 mJ/K'^mol. in very good agreement with 
Ref. [18]. Thus, the additional mass renormalization due 
to the spin-fluctuations exchange is fundamental to rec- 
oncile ARPES and specific-heat measurements. We note 
that the effective masses obtained until now in the 1111 
family by means of de Haas- van Alphen experiments are 
significantly smaller than the values reported in Table 
I.''^ This can be attributed to smaller coupling strength, 
consistently with the smaller Tc values of the fill pnic- 
tides (see also Ref. 44). Recently, de Haas van Alphen 
experiments in BaFe2Asi_j;Pa; show a tendency of con- 
siderably increase of the mass enhancement for this 122 
compound close to the optimal Tc ~ 30 K value. Fur- 
ther de Haas- van Alphen experiments are required to es- 
tablish the correlation between the mass enhancement 
and the values that our analysis suggests. 

To complete the analysis of the specific heat, we com- 
pute numerically the free energy difference Ai^i per band 
between the superconducting and the normal state ac- 
cording to the expression given in Ref. [39], namely 



AF,{T) 



Zfin) 



y^c.2+Af(n)-K| 



(29) 



where Zf, Zf are the Z-renormalization functions for 
the i-band calculated in the superconducting and in the 
normal state respectively, and we evaluate the specific- 
heat difference as ACy,, = -T /^Fi / dT^ . In Fig. 5 
we report the temperature dependence of the ratio of the 
specific heat to temperature for each band, along with the 
total one 75 -7Ar = ACv/T = {Y,i ACv,;)/T. It should 
be stressed that, given the interband nature of the scat- 
tering, the decomposition of the total specific heat ACy 
in single band contributions AC'v,i is purely formal since 
all bands are coupled and the free-energy of one band de- 
pends implicitly on the properties of all the other ones. 
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In agreement with what expected in multiband supercon- 
ductors (see for example Ref. [48,49]), at low temperature 
the increase of 75 is controlled by the quasiparticle exci- 
tations across the smaller gap, i.e. A^g in our case, while 
at higher temperature also quasiparticles excitations in 
the bands with larger gap are thermally activated. The 
presence of multiple energy scales for the quasiparticle 
excitations can lead to the presence of humps in the in- 
termediate temperature range between T = and Tc, as 
observed for example in MgB2 compounds. ^^'^^ In our 
case, such hump is most evident in the contribution of 
the band /3 which is related to the smaller gap. Such 
hump is not however clearly reflected in the total specific 
heat since it is partially compensated by the depletion 
in the other band contributions. The physical origin of 
such depletions is however questionable since it would 
give rise to a negative contributions to the specific heat 
for some band, and it has been argued that it is related 
to the neglecting of superconducting effects in the boson 
propagator itself, as a consequence of the electron-boson 
renormalization.^° The restoring of a physical positive 
contribution of the specific heat for all the bands could 
then make the hump of the /3 band evident also in the 
total specific heat. It is worth noting that it is not clear 
yet if such a hump is visible in the experiments, since 
it is observed in Ref. [17] but not in Ref. [18], where 
a temperature profile remarkably similar to the calcu- 
lations shown in Fig. 5 is reported. Thus, more the- 
oretical and experimental work is required to establish 
the real temperature profile of the specific heat in pnic- 
tides. Finally, we estimate the jump of the specific heat 
at the transition as ^Cy jTc — 72.5 mJ/K^mol, so that 
ACv HTc = 1.45. Both are in good agreement with the 
experimental estimates of ACy /Tc = 98 mJ/K^mol and 
ACv/jTc = 1.5.^''' As far as ACy/^Tc is concerned, it 
must be noticed that even though this estimate is ap- 
parently near to the single-band BCS value 1.43, actu- 
ally in the BCS multiband case the ratio ACy llTc is no 
more universal, so that the experimental result cannot 
be taken as indicative that pnictides are in the weak- 
coupling regime. 

The effect of the multiple gaps is present also in the 
temperature dependence of the superfluid density n^, 
which is experimentally accessible from measurements of 
the penetration depth Al, through the relation A^^ = 
ATTe^ns/mc^ . In two dimensional systems one defines 
conventionally an energy scale Js associated to the su- 
perfluid density n^'* = Ugd of each plane such that: 



^^-^-wAi^^^[^J-^«-^^AiIai^ 

(30) 

where d = 6.6 A is the interlayer spacing. Within the 
Eliashberg approach and in the clean case, can be 
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FIG. 6: (Color online) Temperature dependence of the super- 
fluid density in the four bands within Eliashberg theory. 

computed in each band as:'^^ 

(31) 

Since we are using a parabolic approximation, at T = 
the superfluid density coincides with the carrier density 
rii in each band, so that Eq. (31) reduces to the standard 
formula 

J:(T = 0) = ^ = ^. (32) 

For the band parameters listed in Tables I-II, we can then 
estimate at T = = J2i ^1 — 301 K, while the un- 
renormalized value (obtained with the ARPES band val- 
ues but without taking into account the spin-fluctuations 
induced mass enhancement) would be Js — 700 K. 

The temperature dependence of Jg is reported in Fig. 
6, and it does not differ considerably from the profile ob- 
tained within a simpler multiband BCS approach. ^-"^'^^ In 
all the bands the J|(T) have a flat temperature depen- 
dence at low T, typical of exponentially activated quasi- 
particle excitations across the constant superconducting 
gaps. In the /3 band the deviations from the single-band 
case are more pronounced, due to the low A^/Tc value. 
This anomaly is reflected also in the total superfluid re- 
sponse, that differs from the standard single-band case. 
We notice that the a and 7 bands have almost the same 
normalized profile Js{T)/ Js{0), due to the fact that the 
gaps have approximately the same value in these two 
bands. Thus, the curve in Fig. 6 does not differ qual- 
itatively by a BCS two-band calculation, performed as- 
suming that the band with the large gap contributes to 
63% of the total superfluid density. This is the reason 
why two-band BCS fits (implemented with non-BCS val- 
ues of Ai/Tc) work usually quite well in the comparison 
with the experimental data.^^''^^ 

The comparison of our predictions with the experi- 
mental data is a quite delicate issue due to the pres- 
ence of many extrinsic effects which can spoil a robust 
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Ref. 


a:: 


Tc (K) Xl 


(0) (/im) 


Js(0) (K) temp. dep. 


21 


0.4 


38 


0.231 


307 


exp 


22 


0.45 


30 


0.569 


51 


exp 


23 


0.5 


37 


0.298 


184 


pow 


24 


0.4 


37 


0.208 


380 




25 




32 


0.327 


153 


exp 


26 


0.45 


30 


0.600 


45 


pow 



TABLE III: Summary of some superfluid density measure- 
ments in BKFA compounds in the literature. The last column 
indicate the best fit of the temperature dependence at low T, 
where 'exp' stays for (two-band) exponential fit, and 'pow' 
stays for power-law. 

experimental determination of JsiT). On one hand, 
indeed, measurements of microwave surface-impedance 
have direct access only to \l{T) — Ai(O),^'''^^ so that 
the determination of the A^^(T) profile usually requires 
the separate knowledge of Ai(0), even if in some cases 
one can directly access the normalized penetration depth 
A|(0)/A|(T).2O On the other hand, the /iSR measure- 
ments are also quite delicate because the signal due to 
the screening supercurrents must be disentangled from 
the signal due to magnetic domains, an issue particularly 
delicate in those samples where superconductivity coex- 
ists with residual magnetic ordering. A third problem 
comes from the presence of disorder, which is particu- 
larly severe in pnictides due to the s± symmetry of the 
order parameter, so that inter-band impurity scattering 
is pair-breaking and acts in the same way as magnetic 
impurities in a conventional single-band s-wave super- 
conductor. As a consequence, disorder can lead to a 
change of the low-temperature superfluid-density profile 
from the exponential behavior to a power-law behavior in 
the strong-impurity limit. ^^'^'^ When this is the case, one 
observes also a strong suppression of Js(0) with respect 
to the clean-limit estimate based on Eq. (32).^^ To give 
a hint about how much spread of the data is present in 
the literature we summarize some results for hole-doped 
BKFA in Tab. III. Note that the exact doping of the 
samples is not always available, and that for the same 
nominal doping can be sensibly different, due to the 
different level of disorder. We also included an estimate 
of \l (0) done in Ref. [24] from optical-conductivity data. 
To make a comparison with our clean-limit estimate (32) 
of Js (0) we should then disregard the data characterized 
by a power-law behavior at low temperature, signature 
of strong interband impurity scattering, and the data 
taken for samples with considerably different Tc- The 
best candidates are then the data from Refs. [21] and 
[24], which are taken in samples with the same doping 
level X = 0.4 and same Tc = 37-38 K. These measure- 
ment give Js(0) = 307 — 380 K, which is in very good 
agreement with our Eliashberg calculations. 

Finally, we would like to mention that despite the 
spread of experimental data, in pnictides the T = value 



of the superfluid density is approximately consistent with 
the value of the Fermi energy, apart the not too large 
renormalization effects discussed here. As it was em- 
phasized already in Ref. 11, this implies generically low 
values of the superfluid density in pnictides, simply due 
to the low value of the density of carriers in each band, 
see Eq. (32). This result must be contrasted with the 
case of cuprates, where the density of electrons n is large 
(of order of n ^ 1 — x, where x is the number of doped 
holes/electrons), but nonetheless the superfluid density 
rig is small, and scales approximately with x, so that 
Js ^ xt, where i is a typical hopping energy scale. '''^ For 
this reason, the validity of the well-known Uemura plot,^^ 
i.e. the scaling of with instead of the gap value A 
in underdoped cuprates, does not signal any analogy be- 
tween the two classes of materials. Indeed, in pnictides 
Tc cx Jg a A, with a small Jg simply due to the fact 
that Fermi energy is small, in cuprates the suppression 
of Js with respect to A is due to the proximity to the 
Mott insulator, and the scaling of Tc with Jg can suggest 
a predominant role of phase fluctuations.^'^ 



V. CONCLUSIONS 

In the present work we propose an intermediate- 
coupling Eliashberg multiband approach as an appropri- 
ate description of low-energy properties of pnictides. As 
a starting point we use the bands measured by ARPES, 
where an overall factor two of renormalization of the 
bands with respect to LDA is observed, ^'^^"^^ which orig- 
inate from correlations and cannot be described by the 
coupling to a low-energy bosonic mode (wq — 20 meV) 
discussed here within the Eliashberg theory. We focus in 
particular on BKFA systems, where the multiband struc- 
ture is accompanied by a pronounced anisotropy of the 
Fermi-pocket sizes of the hole bands, with an inner hole 
pocket almost nested to the electron one trough the an- 
tiferromagnetic Q vector of spin fluctuations. We have 
shown that the simultaneous observation of two similar 
gap values in these bands suggests that the predomi- 
nant pairing is an interband one, as mediated by the 
spin fluctuations between the set of hole pockets and the 
set of electron pockets. By comparing the calculations of 
the gaps with the experimental data we have estimated 
the value of the interband coupling, and we calculated 
the corresponding low-energy renormalization in several 
spectral and thermodynamical properties. In particular, 
we showed that a single set of parameter values can ex- 
plain in a consistent way the data on the specific heat^'' 
and on the superfluid density,^^'^** and we can predict the 
appearance of low-energy kinks in the band dispersions, 
that are not always resolved in the experiments. Consis- 
tently with the mean-field character of Eliashberg theory 
we overestimate the critical temperature, which in real 
systems is reduced by superconducting fluctuations,^'' 
whose relevance has been recently pointed out in the con- 
text of paraconductivity measurements."^^ Our analysis 
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questions the possibility of an extreme strong-coupling 
estimate as the one proposed recently in Ref. [36], be- 
cause in this case one would find a huge mass renormal- 
ization at low energy, that is in disagreement with the 
experimental measurements of various thermodynamic 
quantities. Moreover, we have clarified that one must 
resort to a full numerical Eliashberg calculation, an issue 
that has been overlooked in previous studies of multiband 
models with dominant intcrband interactions.^^ Indeed 
we have explicitly shown that in this case the McMillan- 
Allen-Dynes-like approximate expansion^^ fails already 
at coupling values A > 0.2 — 0.3, well below the ones 
relevant for pnictides. 

While the absolute values of the spectral and thermo- 
dynamic properties can be captured only within a four- 
band Eliashberg theory, we have shown that the tem- 
perature dependence of the same quantities do not show 
significant deviations with respect to a two-bands BCS- 
like behavior, once the renormalizcd parameters are used. 
This explains the success of many two-band BCS fits 
proposed in the literature to reproduce the experimen- 
tal data. However, these fits must not be taken as in- 
dicative of the success of a two-band BCS theory, that 
would instead completely fails both from the qualitative 
and quantitative point of view in explaining the physics 
of pnictides. 

While in the present work we focused on BKFA com- 
pounds, our results can be extended to other families of 
pnictides, once that the proper modifications due to the 
different nesting properties of the various Fermi pock- 
ets are taken into account. An interesting example is 
provided by recent ARPES reports in electron-doped 
BaFei.85Coo.i5As2 (Tc = 25.5 K), where around the 
r point only the /3 band crosses the Fermi level, and 
= 6.6 and A-y = 5 meV.^^ In this case, the almost 



nested bands are the /3 pocket and the electron pock- 
ets However, if the bands have the same DOS mea- 
sured in hole-doped compounds and reported in Table 
I, the considerable DOS anisotropy between these bands 
can explain why, even in the presence of the moderately 
large interband pairing suggested by the Ai/Tc values, 
the gaps still differ by 30%. Our multiband Eliashberg 
scheme, with nesting-modulated pairing strength, seems 
a suitable approach to be used to explain the matcrial- 
and doping-dependent hierarchy of the gaps in these pnic- 
tides as well. However, more experimental information 
on the high-energy band renormalization would be re- 
quired to get more quantitative results. Indeed, also 
in BaFei.85Coo.i5As2 one observes the existence of a 
mass renormalization beyond LDA at energy scales much 
higher that the one where spin fluctuations are active. 
A quantitative estimate of these effects would help com- 
paring the overall mass renormalization with recent mea- 
surements of the speciflc-heat in Co-doped BFA,^^ where 
7 seems to be smaller than what found in K-doped crys- 
tals. Indeed, as we suggest in the present work, such 
a comparison is crucial to elucidate the dichotomy be- 
tween high-energy and low-energy mass renormalization 
effects. Thus, further theoretical and experimental inves- 
tigation in this direction can certainly help understanding 
the physics of superconducting pnictides. 
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